PROCEEDINGS 

OF SCIENCE 



o 
o 

(N 

Or 
<D 

r- 



Oh! 

(N 
> 
m 

(N 

^t- 

ON 
O 

o 



Hamiltonian Lattice QCD near the Light Cone 



H.J. Pirner*" c , D. Griinewald" ,E.-M. Ilgenfritz* and E.V. Prokhvatilov" 

" Institut filr Theoretische Physik, Universitdt Heidelberg, Germany 
b Institut fur Physik, Humboldt-Universitdt zu Berlin, Germany 
1 Max-Planck-lnstitut fiir Kernphysik Heidelberg, Germany 
d Department of Theoretical Physics, St. Petersburg University, Russia 
E-mail: |pirner@tphys . uni-heidelberg . de 



d . qruenewaldtjjtphys . uni-heidelberq . de 



ilgenf ri@physik . hu-berlin . de 



Evgeni . Prokhvat@pobox . spbu . ru 
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by standard methods. For this Hamiltonian we determine variationally an approximate ground 
state wave functional in the light cone limit. 
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1. Introduction 

Although the Hamiltonian is not Lorentz invariant, the light cone Hamiltonian [jj], ||] offers 
the advantage of being boost invariant and has - naively interpreted - a trivial vacuum. On the 
other hand, one would be surprised if QCD looses its nonperturbative vacuum structure in the 
light cone limit. Probably much of the complicated vacuum structure of QCD is hidden in the 
constraint equations appearing in light cone QCD. Remarkable progress has been made in light 
cone QCD with a color dielectric lattice theory as a starting point [^, This approach is based 
on "fat" links which arise from averaging gluon configurations [Q]. With this method the spectrum 
of glueballs and the pion light cone wave function have been calculated [Q]. On the light cone one 
is prevented from approaching the continuum limit, since the effective Hamiltonian for the link 
matrices M G GL(N) approaching U G SU (N) is not known. This is the reason why we propose to 
formulate QCD near the light cone using SU (N) link variables. Near light cone time plays a similar 
role as ordinary Minkowski time, therefore we can follow the conventional method of the transfer 
matrix. The transversal field strengths are increased in magnitude due to the boost into the vicinity 
of the light cone whereas the longitudinal fields remain unchanged. Constraint equations arise in the 
light cone Hamiltonian framework which enforce the "equality" of the transverse chromo-electric 
and chromo-magnetic fields = F" k . The lattice Hamiltonian density depends on an effective 
constant which represents the product rj = £, r\ of the asymmetry parameter t, = a_ /aj_ and the near 
light cone parameter r\ . If one chooses r\ = 1 and lets t, — » one obtains a deformed system which 
is squeezed in the spatial a_ direction, if one uses £ = 1 and lets T7 — > one obtains the light cone 
limit. This equivalence has been found by Verlinde and Verlinde ^ and Arefeva These authors 
have proposed to implement such asymmetric lattices in order to study high energy scattering. This 
has motivated us to proceed further in this way. In the work of Balitsky ^ one approaches the light 
cone from time-like distances which is close to the scattering experiments. However, we approach 
the light cone from space-like distances. The asymmetric lattice Hamiltonian itself is not usable for 
Monte Carlo methods since the electric field strengths i.e. the momenta appear linearly. Because 
of the translational invariance of the vacuum we can add a term 1/t] 2 P_ to cancel the unwanted 
terms. Naively this amounts to returning to an effective lattice Hamiltonian which is proportional 
to the energy in ordinary Minkowski coordinates. This is a reasonable procedure to search the 
groundstate of the vacuum. Applications of the light cone coordinates in finite temperature field 
theory have followed the same route [Q]. 



2. The QCD Hamiltonian near the light cone 



We introduce near light cone (nlc) coordinates, first proposed by [jlfj 
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The transversal coordinates x 1 and x 2 remain unchanged. The near light cone parameter v\ may 
be interpreted as parameterizing a Lorentz boost into a frame which is moving with velocity j8 = 
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(1 — T] 2 /2) /(l + T] 2 /2) along the longitudinal direction relative to the laboratory frame. Then, the 
nlc energy p + and longitudinal momentum p_ expressed in terms of the laboratory energy Ei a b and 
longitudinal momentum p? , are given by 
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(2.2) 



The second relation in eqn. (2J2) shows that large longitudinal momenta p? , become accessible by 
a nlc lattice with a cut-off l/a_. 

The definition of nlc coordinates eqn. (2T) induces the following metric: 

/0 1 \ /t] 2 l\ 
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1. This defines the scalar product 
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Note, that the metric has off-diagonal terms which implies that there are terms mixing temporal 
and longitudinal spatial coordinates in the scalar product. This yields a form of the pure gluonic 
Lagrange density eqn. ( |2.5[ ) which has severe consequences for a numerical treatment 

JSf = ' 



(2.5) 



The Lagrange density is linear in one of the temporal field strengths, namely F" k F" k . Therefore, 
a standard Monte Carlo sampling of the Euclidean path integral does not work for nlc coordinates 
and we rather use a Hamiltonian formulation. 

The energy momentum tensor in its most general form is given by 

- 8^ 



It defines the Hamiltonian density Jif 
as 
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T\ and the longitudinal momentum density 
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This form of the local integrand for the generator of longitudinal translations is not manifestly 
gauge invariant. However, if one uses Gauss' law and the definition of the field strengths tensors 
one can rewrite in a symmetrized form 



(2.8) 
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In order to solve the Hamiltonian we are interested in translation-invariant states which are eigen- 
states of the longitudinal momentum operator, i.e. with eigenvalue equal zero. In vacuum, with 
light cone momentum P = 0, we can add (1/t] 2 ) P to define an effective Hamiltonian density 
M'eff which is only quadratic in momenta: 
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In a forthcoming paper we show how to derive [ ]11[ ] the effective lattice Hamiltonian with the 
coupling constant X =4/g 4 using the transfer matrix: 
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For f} = 1 this effective Hamiltonian is very similar to the traditional Hamiltonian used in equal 
time lattice theory. It differs only in the potential energy terms for the plaquettes. Instead of 
the usual Tr[l — Re(f/_^)] term resembling the field strength squared in the naive continuum limit, 
the nlc Hamiltonian has the form (Tr[a"/2 lm(U-k)]) 2 which corresponds to the plaquette in the 
adjoint representation and which yields an additional Z(2) symmetry in comparison with the full 
lattice Hamiltonian [11]. The light cone limit 77 — ^ enhances the importance of transverse elec- 
tric and magnetic fields without generating unwanted linear terms in the momenta. The resulting 
vacuum solution should be a plausible extrapolation of the vacuum solution of QCD. 



3. Variational optimization of the ground state wavefunctional 



We analytically solve the effective lattice Hamiltonian for the ground state wavefunctional 



both in the strong and weak coupling limit [11]. Both solutions can be described by a product 
of single plaquette wavefunctionals for f) sufficiently close to one. In order to cover the whole 
coupling range, we make a variational Ansatz of the ground state wavefunctional which is given 
by a product of single plaquette wavefunctionals with two variational parameters p and 5 and with 
the normalization constant N 
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With this normalized wavefunction we optimize the energy expectation value £b(p, §) of the effec- 
tive Hamiltonian for fixed values of the coupling A and the near light cone parameter tj . 

In fig. ([I]) and fig. (Q), we present the variationally optimized wavefunctional parameters po 
and 8q rescaled by a factor 1 / vX such that they approach a constant in the assymptotic weak 
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Figure 1: Optimal wavefunctional parameter po(X,fj) as a function of A obtained from the simulation 



on a Nf x N- = I6 i lattice for three different values of f) 2 . The red shaded area corresponds to the phase 



transition region for all values of f] 2 . The dotted lines show the predicted analytical strong coupling behavior. 
The arrows indicate the expected asymptotic behavior for weak coupling which is proportional to \fk, i.e. a 
constant independent of X in the plot. The solid lines show the actual analytic parameterizations in the weak 
coupling regime (c.f. eqn. (3.3)). 



coupling region, i.e. A — » °°. The statistical uncertainty on the variational parameters is typically 
5%. It becomes larger in the region where the Hamiltonian with the (— &)-plaquette in the adjoint 
representation induces a phase transition related to the Z(2) symmetry. This region is indicated by 
the red shaded area in the figures. In principle only couplings in the weak coupling region above 



A = 7 are physically meaningful where the artificial Z(2) symmetry in eqn. ( |2.10| ) is spontaneously 
broken. 

The discussed analytical strong and weak coupling solutions [JTTJ] yield the following estimates 
for po and 8q in these limits 



Po 
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for A << 1 
for A >> 1 



, /. for A << 1 

A fj 2 7fj(0) for A >> 1 



0.038 

M°) - —~-,forf)^0. (3.2) 

The variationally determined parameters are in good agreement with the analytic predictions in the 
strong coupling regime which are represented by the dotted lines in fig. ([[]) and fig. @. However, 
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Figure 2: Optimal wavefunctional parameter 5o(X,fj) as a function of X obtained from the simulation 
on a Nj_ x N- = 16 3 lattice for three different values of fj 2 . The red shaded area corresponds to the phase 
transition region for all values of fj 2 . The dotted lines show the predicted analytical strong coupling behavior. 
The arrows indicate the expected asymptotic behavior for weak coupling which is proportional to \fX, i.e. a 
constant independent of X in the plot. The solid lines show the actual analytic parameterizations in the weak 
coupling regime (c.f. eqn. (3.3)). 



in the weak coupling regime the optimal parameters differ from their analytical estimates eqn. ( |32| ) 
which are indicated by the arrows in the plots. Both analytic predictions disagree with the optimized 
values stronger for decreasing values of 77 . This is natural, since the light cone limit fj — > builds 



up correlations among plaquettes separated along the longitudinal direction Qll[]. The parameters 
optimizing our product of single plaquette wave functionals can only effectively describe these 
correlations. 

In the physical relevant coupling region beyond X = 7, it is possible to fit the variationally 
optimized wavefunctional parameters to the following parameterization 
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(3.3) 
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A good fit of the parameters coj,...,C5j minimizing x 2 in the range X € [10,95] and fj G [0.15, 1] is 
possible and yields the coefficients tabulated in tab. ([[]). The result of the fitting procedure eqn. (3.3) 
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co,i 


c l,i 


C 2J 


c 3,r 


C4,i 




P 


0.90 


-1.74 


0.72 


8.12 


-0.40 


-0.27 


8 


0.95 


0.93 


-1.21 


-6.44 


-0.83 


0.64 



Table 1: Coefficients of eqn. (3.3) obtained from least square minimization. 

is shown by the solid lines in fig. ([[]) and fig. @. Having a parameterization of the ground state 
wavefunctional in dependence of the nlc parameter fj at hand, we plan to extrapolate to the light 
cone and finally calculate hadronic cross sections by simulating how a color dipole moving near 
the light cone hits a neutral hadron localized at x~ =0 [12]. The color dipole can be represented 
by a longitudinal-transversal Wilson loop elongated in x~~ direction and the simplified target by 
a transverse plaquette. Varying the impact parameter one can sample the correlation function of 
the two gauge-invariant objects and thereby obtain the profile function. A necessary prerequisite 
of such a calculation for different velocities of the dipole is that the lattice constant in transverse 
direction stays constant for different fj values, in order to have a reliable transverse length scale. 
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